Computer Science and Information Technologies 
Vol. 1, No. 3, November 2020, pp. 126~134 
ISSN: 2722-3221, DOI: 10.11591/csit.v1i3.p126-134 im) 126 


An optimized rubber sheet model for normalization phase of 
iris recognition 


Selvamuthukumaran.S!, Ramkumar.T’, Shantharajah SP? 
'Faculty of Computer Applications, AVC College of Engineering, India 
3S$chool of Information Technology and Engineering, VIT University, India 


Article Info ABSTRACT 

Article history: Iris recognition is a promising biometric authentication approach and it is a 
. very active topic in both research and realistic applications because the pattern 

Received Jan 23, 2020 of the human iris differs from person to person, even between twins. In this 

Revised May 27, 2020 paper, an optimized iris normalization method for the conversion of segmented 

Accepted Jun 16, 2020 image into normalized form has been proposed. The existing methods are 


converting the Cartesian coordinates of the segmented image into polar 
coordinates. To get more accuracy, the proposed method is using an optimized 
Keywords: rubber sheet model which converts the polar coordinates into spherical 
coordinates followed by localized histogram equalization. The experimental 
result shows the proposed method scores an encouraging performance with 
respect to accuracy. 


Biometrics 
Histogram equalization 
Iris normalization 


Iris recognition 
Spherical coordinates This is an open access article under the CC BY-SA license. 


Corresponding Author: 


Selvamuthukumaran.§S, 

Faculty of Computer Applications, 
AVC College of Engineering, India. 
Email: smksmk@gmail.com 


1. INTRODUCTION 

The progressive digital society demands secure identification techniques for the development of 
biometric system in diverse fields [1]. The iris biometric systems are becoming commonly implemented as one 
of the best ways to certainly identify people [2]. The iris is an outwardly and colourful organ close the pupil of 
the eye [3]. Figure 1 shown the structure of the eye, which demonstrates the exact position of the iris and its 
surrounding things. The property of the iris guarantees that even equal twins have uncorrelated iris details. 
Thus, the exclusivity of every iris, including the couple possessed by single individual, parallels the exclusivity 
of every fingerprint regardless of whether there is a common genome [4, 5]. The iris involves of several 
irregular small blocks similar to freckles, stripes, furrows, coronas, etc. In addition, the texture divisions in the 
iris are arbitrary. The separate merits of the iris cause its high consistency for individual identification [6]. 
Hence the method of iris identification becomes an exciting research point in current years [7]. The first step 
in iris recognition is the process of capturing the image of an eye [8]. During stage of the image acquisition, a 
camera is used to record iris textures. The captured image is further processed for the localization stage, where 
the location of the iris is carried out, followed by the segmentation phase. In segmentation, irises are separated 
from the eye [6]. Next to this, a normalization of the iris image is done using various noise removal techniques, 
and the images are stored in the binary format [9]. The process of the iris recognition system continues with 
the feature extraction stage, where the extractions of various features of irises are identified [10]. In verification 
process, matching will be done where the acquired image is compared with the images stored in the database 
[11]. A typical iris recognition system which includes above stages is shown in the Figure 2. 
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Figure 2. Iris recognition system 

In the phase of normalization, segmented iris image has been prepared for the feature extraction 
process. Though, several of algorithms [1, 5, 11] are available for the normalization phase of the iris 
recognition, Daugman [12] and Wildes [13] algorithms are significant. To compensate the differences in the 
magnitude of the pupil, Daugman [12] proposed a rubber sheet model for normalization. In this way, some of 
the irises still may be sealed by eyelids or eyelashes, even when the inner and outer boundaries of the iris are 
reached. Wildes’ system is a patent of the iris recognition system which employes the gradient-based Hough 
transform to choose the two circular boundaries of an iris. It has big computational charge, since it examinates 
among all of the possible candidates. Besides, the algorithm’s accuracy significantly deteriorates while dealing 
with the noisy data. The above two normalization techniques transform Cartesian coordinates in to polar 
coordinates for un-wrapping the iris texture into a fixed size rectangular block. In polar coordinates, the 
distance and angular position have extremely affected iris images with respect to the camera. Additionally, 
illumination has a direct influence on the pupil size and sources non-linear distinctions of the iris patterns [14]. 

Moreover, the pupil boundary and limbus boundary are usually two non-concentric contours. The 
non-concentric condition leads to different choices of reference points for transforming an iris into polar 
coordinates, which is not suitable for iris images with more noise. Hence, it is necessary to adopt an 
optimization technique for transforming the iris image to compensate these variations. In this paper, we have 
proposed an optimized normalization technique based on rubber sheet model which converts Cartesian 
coordinates into spherical coordinates instead of polar coordinates and perform localized histogram 
equalization for further accuracy. The rest of the paper has been organized as follows; in section 2, the 
improvised model for iris normalization has been illustrated with the inclusive steps. The experimental 
investigations of the proposed work with the help of CASIA database [15] is explained in section 3, where 
concluding remarks and feature research directions are summarized in section 4. 


2. AN OPTIMIZED RUBBER SHEET MODEL FOR IRIS NORMALIZATION 

The normalization stage aims to gain invariance to size, position and pupil dilatation in the segmented 
iris region. Most of the methods use Daugman [12] rubber sheet model considers the possibility of pupil dilation 
and the appearance of different sizes in different images. For this purpose, the coordinate system is altered by 
un-wrapping the iris and plotting all the points within the edge of the iris into their polar coordinates. The 
Cartesian to polar transformation is defined for the continuous form of images. However, in the discrete form, 


An optimized rubber sheet model for normalization phase of iris recognition... (Selvamuthukumaran.S) 


128 0 ISSN: 2722-3221 


the transformation encounters problems in which the polar samples do not entirely match the Cartesian samples 
which results in excessive interpolation and in some cases loss of information would be the result. The proposed 
method is the improvised method of Daugman’s rubber sheet model and it improves the accuracy by advocating 
a two-step process. At First, the coordinate system of the iris region is converted from Cartesian into the 
spherical. Then, the remapping process is done for the normalization. Localized histogram equalization is 
applied on the iris region for further normalization. The various inclusive stages of the proposed work are 
illustrated in the Figure 3. 
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Figure 3. Flowchart of proposed work 


2.1. Conversion of cartesian coordinates into spherical coordinates 

In the spherical coordinate system, the position of a point is specified by three coordinates; (r, 8, p) 
where ‘r’ is the radial distance from a fixed origin, ‘0’ is the elevation angle of that point from a plane, and 
‘gp’ is the azimuth angle of its orthogonal projection on that plane from a fixed direction. The elevation angle 
is often replaced by the inclination angle measured from the azimuth; the direction perpendicular to the 
reference plane. The radial distance is also named the radius or radial coordinate and the inclination is named 
as colatitude. The coordinate conversion is represented in Figure 4. During co-ordinate conversion, the centre 
of the pupil is considered as location point and the radial vector circle has been considered in the iris region. 
The coordinate conversion is performed to obtain the invariance of iris size, position and different degrees of 
pupil dilation. It produces the iris regions, which have the same constant dimensions, with the intention that 
two images of the same iris should have characteristic features under different conditions at the same spatial 
place. The proposed technique is useful in bringing iris images into a standard fixed resolution, which simplifies 
the feature extraction process. In the adjusted rubber sheet model, the centre of the pupil is used as the reference 
point; and the radial and azimuth vectors cross the iris region. The spherical coordinates, (r, 0,@) of a point 
can be obtained from its Cartesian coordinates, (x, y, z) by the following equations. 
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Figure 4. Coordinate conversion 


r= Vx*+y?+z? (1) 


g =atan (y,x) (2) 
6 =arccos —_ (3) 


Where atan (y, x) is a variant of the arctangent function that returns the angle from the x-axis to the 
vector (x,y) in the full range(—z,7). The formulas assume that the two systems have the origin, then the 
spherical reference plane is the Cartesian 'x — y’ plan ‘@’ is the inclination from the ‘z’ direction and the 
azimuth angles are measured using the coordinates. After converting the data points of the iris region into 
spherical coordinates, remapping has been done to rescale the points on the angle around the inner and outer 
boundaries. The remapping of the iris region from ‘(x,y,z)’, the Cartesian coordinates to the normalised 
spherical representation can be modelled as in the following equation. 


(xr, 8,9), Y7, 8,9), 2(7, 8, P)) =1(7, 8, 9) (4) 


Where 'I(x,y,z)’ is the iris region image, ‘(x,y,z)’ are the original Cartesian coordinates, and 
‘(r,@,@)' are the corresponding normalized spherical coordinates. It also represents the coordinates of the pupil 
and iris boundaries along the ‘0’ direction and furnished in the following equations. 


In(%,¥,2) = Io(X,y, 2) (5) 
x = X,(8,9) + (Xi, 9)- xp(8, 9))— (6) 
Y =Yp(9,9) + (iO, 9)- Yp(8,.9))— (7) 
z= 2,(8,9) + (2(8,9)- 2p(9, 9))— (8) 


Where 0 = 27 - 
Here ‘I,’ is the ‘'MXN’normalized image and 'x,(6,@),¥p(8,~),Z,(8,~)' and 
'x:(0,~), vi(8,), Z:(8,)' are the coordinates of the inner and outer boundary points. ‘0’ is the inclination 


from the ‘Z’ direction in the original image, ‘Ij’. This process ends with the normalised iris image, which is 
further forwarded for histogram equalization with the view of getting better accuracy. 


2.2. Localized histogram equalization 

The normalized image, resultant of the previous phase may have low intensity contrast and improper 
brightness due to the position of the light sources and other issues. In such situations, to achieve more accurate 
results, histogram equalization has been chosen. This method is a contrast enhancement technique with the 
objective to obtain a new enhanced iris image with a uniform histogram. It distributes the histogram of the 
pixel levels in a uniform manner without any loss. It has been implemented for improving the look of a noisy 
image by amplifies the global contrast of the iris image. The histogram, ‘H[x]’ contains pixels with the value, 
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'x'. The cumulative density function of the histogram ‘cdf [x]’ contains the number of pixels with the value, 
'x' or less is given by the following equation. 


(cdf [x]) = H[0] + H[1] + H[2] + ---+ A[x] (9) 


It should be applied to each pixel in the iris image and replace the existing value with the calculated 
value. The histogram of a digital image in the range of [0,K —1] is a discrete function as given in 
the following equation. 


p(T) = “* (10) 
Where ‘ 7 ’ is the k" gray level, is the number of pixel in the image with that gram level, ‘n’ is the total number 
of pixels in the image, and k = 0,1,2,...,K-1. p(%) Gives the probability of occurrence level, r,. 

Consider the following sample original entire eye shown in the following Figure 5. The contrast of 
the iris image is amplified, mainly when the used data of the iris image is represented by near contrast values. 
The intensity would be better spread on the histogram, through this fine-tuning. It changes the areas of lower 
contrast to rise a higher contrast without distressing the overall precision of the irises. The feature of interest 
in the image may require enrichment locally. Histogram equalization reaches this by efficiently spreading out 
the most frequent intensity values. By using this method, the histogram of the resultant image is as flat as 
possible. The resultant histogram equalized image is shown in the Figure 6. 


Figure 5. Original image before localized histogram equalization 


Figure 6. Image after the localized histogram equalization 


It consists of applying localized histogram equalization independently to capture the iris region in the 
eye image since most small regions are very self-similar. If the image is made up of discrete regions, most 
small regions lie entirely within one or the other region. This method attempts to equalize the number of pixels 
which tend to flatten and raise an image’s histogram. Vertical equalization uses only a single column of pixels 
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into the equalization process, whereas horizontal equalization uses a single row of pixels. It modifies the pixel 
intensities for a better appearance of the irises. The proposed method improves the visual appearance [16] of 
an image. 


3. EXPERIMENTAL STUDY 

To achieve a size-invariant sampling of the valid iris pixel points, we have applied the optimized 
rubber sheet model to map the sampled iris pixels from the Cartesian coordinates into the normalized spherical 
coordinates. To effectively implement our proposal, the iris image database CASIA Database Version 3.0 
(CASIA-IrisV3) released by the Centre for Biometrics and Security Research of National Laboratory of Pattern 
Recognition, Chinese Academy of Sciences has been used [15]. The pupil is smaller in the image; however, 
the normalisation process can able to rescale the iris region, so that it has a constant dimension. In our 
experimental study, the rectangular representation is constructed from 12,000 data points in the iris region. The 
remapping is done using the spherical coordinates. The iris shown in Figure 7 is the segmented image of 
S1046R04. 


re 


Figure 7. Segmented iris image of S1046R04 


After the successful conversion of the Cartesian coordinate into the spherical coordinate system, the 
normalised iris is shown in the following figure. Comparing the results from Figure 8 and Figure 9, it has been 
observed that the distribution is shifted towards the higher values, while the peak at the minimum intensity 
remains after the localized histogram equalisation. It recovered some of the apparently lost contrast in an image, 
by remapping the brightness values in such a way as to equalize, or more evenly distribute its brightness values. 
The histogram of the normalised iris of Figure 8 is shown in Figure 10. The unwrapped flat iris has a low 
contrast. This iris is enhanced by applying localized histogram equalization. The resultant image is shown in 
the following figure. The histogram of the normalised iris of Figure 9 after the localized histogram equalization 
is shown in Figure 11. 


Figure 8. Normalized iris 
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Figure 9. Normalized iris after histogram equalization 
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Figure 10. Histogram of normalized iris after conversion into spherical coordinates 


Final 
Figure 11. Histogram of the normalized iris after applying the localized histogram equalization 
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By comparing the histograms in Figure 10 and Figure 11, it has been observed that the intensity 
variations are equally normalised after the localized histogram equalisation. Also output of the proposed 
method for some of the CASIA images are shown in the Table 1. As a result, noises in the flat regions and ring 
artefacts at the edges are normalised, and now the irises are ready for further processing. 


Table 1. Result of the proposed method for some of the CASIA images 
S.No CASIA Image 
1 $1013L02 


2 $1022L04 


The proposed segmentation method normalized the iris region in 658 images out of 670 images, which 
corresponds to a success rate of 98%. The success rate of the proposed method is compared with those of the 
previous methods Daugman’s Rubber sheet model and Wild’s Image registration. The success rate of the 
different normalization methods is shown in Table 2. Table 2 shows the strength of the proposed method which 
works better than the existing methods since 98% of accuracy for the irises and the graphical representation of 
comparison is also provided in the Figure 12. 


Table 2. Success rate of the normalization methods 


Normalization Method No of irises (out of 670) | Success Rate percentage 
Optimized rubber sheet model using spherical coordinates 658 98 
Rubber sheet model 638 95 
Wild’s Image registration 597 89 
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Figure 12. Results of various normalization methods 


CONCLUSION 
In a progressively digital society, the role played by the iris recognition system is a vital one. The 


trends and traits made in security-based technologies demand various new algorithms for quality images. In 
this paper, an optimization technique for the normalization phase has been implemented which performs the 
conversion of the Cartesian coordinates into spherical coordinates and uses localized histogram equalization. 
The results obtained reveal the exact normalization of irises for the noisy images. As a future work, we put 
forth certain improvements in noise removal methods, iris region resize methods towards the better feature 
extraction of iris images, 
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